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Abstract 

We combine density functional theory (DFT) within the local density ap- 
proximation (LDA), the quasiharmonic approximation (QHA), and a model 
vibrational density of states (VDoS) to calculate elastic moduli and sound 
velocities of 7— (Fe^Mg^^SiO^ (ringwoodite), the most abundant mineral 
of the lower Earth's transition zone (TZ). Comparison with experimental 
values at room-temperature and high pressure or ambient-pressure and high 
temperature shows good agreement with our first-principles findings. Then, 
we investigate the contrasts associated with the (3 — >■ 7(Fe x ,Mgi_.c)2Si04 
transformation at pressures and temperatures relevant to the TZ. This infor- 
mation offers clearly defined reference values to advance the understanding 
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of the nature of the 520 km seismic discontinuity. 

Keywords: first principles, ringwoodite, elasticity, transition zone, 520 km 
discontinuity 



i 1. Introduction 



Wadsleyite (/3-phase) and ringwoodite (7-phase) are the high-pressure 
polymorphs of olivine (a-phase), (Fe x ,Mg 1 „ a ,) 2 Si04. These minerals are the 



main constituents of the Earth's upper mantle (UM) ( 



Rinwood 



19921 ) and transition zone (TZ) ( llrifune and Ringwoodl . Il987| ). Under pres- 



1975 



Putnisl . 



e sure, the transformation from olivine to wadsleyite happen s at ~13.5 GPa 



and f r om wadsleyite to rin gwoodite at ~18 GPa near 1600 K (IKatsura and Ito 



198 



Akaogi et al. 



19891 ) . These transformations are associated with two 



major discontinuities in seismic veloci ties in the Earth's interior at a bout 410 
km and 520 km depth, respectively (IRevenaugh and Jordan! Il99l[ ). While 
the first discontinuity is a well characterized and sharp feature in seismic 
data, the second varies considerably with location. 

Experimental and computational approaches have been used to study 
properties of these minerals at in situ conditions and their relationship with 
seismic discontinuities. Mg-end member and Fe-bearing a-phases have been 
widely investigated at simultaneous high pressure and temperature. Compu- 
tational results and experimental data on elastic properties a nd sound veloci- 



ties seem to be consistent with seismic measurements (see e.g 



2010 



Stackhouse et al 



2010 



Li and Liebermann 



Nunez Valdez et al. 



2007, and references there 



in). On the other hand, even though great efforts have been made to ob- 
tain measurements of elastic properties and wave velocities of Fe-free and 



2 



22 Fe-bearing wadsleyite and ringwoodite under high tempera ture anc 



23 sure using either ultrasonic or Bri 



Zha et al. 



2000 



2006 



1997; 



Liu et al.. 2009 



Isaak et al. 



Weidner et al. 



2007 



louin scattering techniques ( 



Sinogeikin et al 



Mayama et all 12004 : 



1984 



Sinogeikin et al 



1998 



2003 



Li et al 



pros- 



1996 



Li and Liebermann 



Isaak et all I201oI:Il1 2003: Higo et al. 



Jackson et al 



2000 



Mayama et al. 



20051 ). results are still limited. Therefore large extrapolations from room con- 
ditions to conditions of the TZ are often used. 

First principles calculations employing the quasiharmonic approximation 
(QHA), valid up to about two thirds of the melting temperature, or molecular 
dynamics (MD) methods, valid near and above melting temperatures, com- 
plement each other and are used to obtain elastic moduli under high-pressure- 
temperature conditions. Calculations of elastic constants using the QHA, 
though computationally less demanding than MD, still required calculations 
of vibrational density of states (VDoS) for each strained atomic configura 



36 tion at seve ral pressures, that is, about 1000 parallel jobs (IDa Silveira et al. 



2008 



2011 



)• 



In this paper we use the analytical and computational approach by 



Wu and Wentzcovitch 



39 (120 111 ) and tested on 



^-(Fe^Mgi^SiC^ (INunez Valdez et al 



periclase-MgO, a— Mg?S i04, and more recently on 



20121 ) . to calculate bulk (K) and 



shear (G) moduli and sound velocities of the 7— (Fe^Mgi-^SiC^ phase. 
This method uses only static elastic constants and phonon density of states 
for unstrained configurations, therefore reducing the amount of computa- 
tional time and resources by one to two orders of magnitude. We then ad- 
dress contrasts across the (3 — > 7(Fe a; ,Mgi_. ! ;)2Si04 transition near conditions 
of the 520 km seismic discontinuity. 
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2. Methodology 



2.1. Computational Details 



49 C alculations based on Den sity Functional Theory (DFT) ( IHohenberg and Kohn 



1964 



Kohn and Shaml. 



proximation (LDA) (jCeperley and Alder 



1965 ) were performed using the local density ap- 



52 generated by the Vanderbilt method (jVanderbilt 



1980). U 



traso ft pseudopotentials 



19901 ) were used to de- 



53 scribe Fe, Si, and O. A norm- conserving pseudopotential generated by the 



von Car method was used for Mg. Furt 



tentials are given in (INunez Valdez et al. 



r er de tails about these pseudopo- 



201ll ). Equilibrium structures of 



ringwoodite (28 atoms/cell) at arbitrary pressures were 



57 able cell-shape damped m olecular dynamics approach (jWentzcovitchl . 



Wentzcovitch et al 



bund using the vari 



1991 



59 code (IGiannozzi et al. 



1993) as implemented in the quantum-ESPRESSO (QE) 



20091 ) . The plane- wave kinetic energy cutoff used was 



40 Ry and for the charge density 160 Ry. The k-point sampling of the charge 
density was determined on a 2 x 2 x 2 Monkhorst-Pack grid of the Brillouin 
Zone (BZ) shifted by (|||) from the origin. These parameters correspond 
to having interatomic forces smaller than 10 -4 Ry/a.u. and pressure con- 
vergence within 0.5 GPa. Dynamical matrices were obtain ed using density 
functional perturbation theory (DFPT) ( IBaroni et all l200lf ) via QE. At each 
pressure, a dynamical matrix was calculated on a 2 x 2 x 2 q-point mesh for 
one atomic configuration only. In principle about 10 other configurations 
should be used as well, but here we are more interested in frequencies with 
strain and the current approximation seems to be sufficiently accurate. Force 
constants were extracted and interpolated to a 12 x 12 x 12 regular q-point 
mesh to produce VDoS. 
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2.2. High-Temperature-Pressure Elastic Theory 

Exploiting the information about the strain and volume dependence of 
phonon frequencies, we determi ne the th e rmal contribution to the Helmholtz 



75 free energy F within the QHA ([Wallace! 



1972|), that is, 



F (e, V, T) = U st (e, V) + \ ^q>™( e > V) + 



q,m 



exp 



q,m 
fthJr 



k B T 



(1) 



where q is the phonon wave vector, m is the normal mode index, T is tem- 
perature, U st is the static internal energy at equilibrium volume V under 
isotropic pressure P and infinitesimal strain e, h~ and k B are Planck and 
Boltzmann constants, respectively. Isothermal elastic constants are given by 

'd 2 G(P,T) 



^ijkl 



deije k i 



with G = F + PV, the Gibbs energy, and i, j, k, I — 1, 
adiabatic elastic constants, one uses the relationship: 

T dS dS 



^ijkl 



(2) 

3. To convert to 
(3) 



VC V deij de H 

where Cy is heat capacity at constant volume, and S is entropy. For or- 
thorhombic crystals, the non-shear elastic constants of Eq. ([2]) are: 

~d 2 F(e,V,T) 



C 



;i-<yp(v,T), 



C?t 17 (V) + C!*(V,T) 



(4) 



while the shear elastic constants are: 



cL, = c^(y) + cf^(v,T). 



mo 



(5) 



Elastic constants Cf^jj and Cf^ can be expressed as functions of the volume 
Griineisen parameters, 7 qim = — <9(lnu; qim )/<9(m V): 



dou, 



q,m 



dV 



(6) 



or the generalization to strain Griineisen parameters: 

due 



(7) 



q,m 



We have used the Wu-Wentzcovitch method (iWu and Wentzcovitchl . l201l[ ) 

89 to compute the thermal contribution to the elastic constants, Cf^j(V, T). 

90 This method allows the computation of thermal elastic constants without 

91 performing phonon calculations for strained configurations with the approx- 

92 imation that strain and mode Griineisen parameters have isotropic distribu- 

93 tion, which is equivalent t o assuming th a t ther mal pressure is isotropic. This 



20071 ) implicit in the QHA calcula- 



is a good approximation (jCarrier et al. 
tion of thermal pressures. 

After obtaining VDoS at several volumes by first principles, average 
strain Griineisen parameters were computed at such volumes and interpo- 
lated in a fine volume-temperature grid, that was then inverted to a pressure- 
temperature grid of 0.1 GPa and 10 K spacing s. Static elastic constants 



previously computed ( jNunez Valdez et al.l . 1201 ll ) were also used. 



Voigt and Reuss bounds of bulk and shear moduli of orthorhombic crys- 
tals at high temperature were calculated using adiabatic elastic constants as 



6 



103 (W att et al. 



1976 



Watt 



19791 ): 



Ky = - [Cn + C22 + C33 + 2 (C12 + C13 + C23)] , 
y 

GV = — [Cn + C22 + C33 + 3 (C44 + C55 + C*66) 
15 

— (Cl2 + Cl3 + C23)] j 

Kr = D [Cu (C22 + C 33 - 2C23) + C22 (C33 - 2Ci 3 ) 
— 2C33C12 + C12 (2C23 — C12) 
+C13 (2C12 — C13) + C23 (2C13 — C23)] 1 1 

G B = 15 {4 [Cn (C 22 + C33 + C 23 ) + C 22 (C 33 + C13) 

+C33C12 — C12 (C12 + C23) — C13 (C13 + C12) 

-C 23 (C 23 + C13)] /-D 

+3(1/C 4 4 + 1/C 5 5 + 1/C66)}" 1 , 
D = Ci3 (C12C23 — C13C22) + C23 (C12C13 — C23C11} 
+C33 (CllC22 — C12) ■ 

104 With Voight-Reuss-Hill averages of elastic moduli, i.e., 



K 



K v + K 



R 



and G 



Gy + G 



R 



105 isotropic sound velocities are given by: 
V P = 



[ K+\G 



Vs 



P 




(8) 



0) 



(10) 



(11) 
(12) 

(13) 
(14) 



106 where p is density, and Vp, Vs, and V$ are compressional, shear and bulk 

107 velocities, respectively. 



108 

109 
110 



3. Results 

We present first-principles results of aggregate properties of Fe-bearing 
ringwoodite at pressures and temperatures relevant to the TZ. All the ap- 



Y-Mg _SiO z 



r(Fe x ,M gi . x ) 2 Si04 




d> 
> 



O 300 K, Li 2003 

□ 300 K, Higo et al. 2006 

X 300 K, Weidner et al. 1 984 

■0,1500,2000,3000 K, Li et al. 2006 



A 300 K, Sinogeikin et al. 2003, x=0.09 
300 K, Higo et al. 2006, x=0.2 



5 10 15 20 25 
Pressure (GPa) 



5 10 15 20 25 
Pressure (GPa) 



30 



Figure 1: (Color online) Pressure and temperature dependence of bulk modulus (K), shear 
modulus (G), compressional velocity (Vp), shear velocity (Vs) and bulk velocity (V$) for 
Fe-free ringwoodite (a,c) and Fe-bearing ringwoodite (b,d). First principles calculations 
within LDA (solid lines) are compared to available experimental data (symbols). Note, 
however, that low-pressure-high temperature calculated trends (dash lines) are outside the 
validity of the QHA. 
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proximations described in the previous section provided an excellent de- 
scription of bulk and shear moduli, and sound velociti es within the v alid 



113 regime of the QHA established for the Fe-free quantities (jYu et al. 



2008|) for 



7-(Fe I ,Mg 1 _ x .) 2 Si04 with x = and x = 0.125 (see Figs. [Hand [2]). 

In the case of the Mg-end member ringwoodite, K increases more rapidly 
than G, as a function of pressure, and decreases faster than G with increas- 



n7 ing t e mperature. At 3 K the agreement b etween experimental data (ILi . 



2003 



Higo et al. 



2006 



Weidner et al. 



34) 



19841 ) and our DFT-results is truly 



excellent for elastic moduli (Fig. [T^l) , an d sound velocit ies (Fig. [Tb). Re 



Li et al 



(120061 ) also compare well 



suits from a molecular dynamics study by 
with our results for K, Vp, and V$ within the QHA limits. Our predicted G 
and Vcj are larg er and smaller, respectively, than molecular dynamics values 



Li et al. 



(120061 ). This difference is primarily caused by the use of the GGA 
approximation in the MD simulation. Nevertheless the general agreement 
is good and it is the only other source to compare aggregate properties of 
7— Mg 2 Si04 at high pressures and temperatures. Results of Fe-bearing ring- 
woodite as function of pr e ssure are shown in Figs. [Tb and[TH. Experimental 



data by 



Sinogeikin et al. 



(120031 ) at room temperature and x = 0.09 are in 

129 excellent agreement with our 300 K curve, though after 10 GPa our K tends 

130 to be larger, while G tends to be smaller (Fig. [Tb). Fe-bearing compres- 

131 sional, shear and bulk velocities are smaller than their Fe-free counterparts 



132 (Fig. [TU). Predictions for x = 0.125 fal 



133 with iron concent rations of x = 0.09 (ISinogeikin et al 



134 (IHigo et al 



in between two experim ental reports 



20031 ) and x = 0.2 



20061 ) . From Fig. [TH one can see that Vp is the most affected by 



135 iron concentration and temperature, while V$ seems to be the least affected 
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Figure 2: (Color online) Temperature dependence of bulk modulus (if), shear modulus 
(G), compressional velocity (Vp), shear velocity (Vs) and bulk velocity (V$) for Mg-end 
member ringwoodite (a,c) and Fe-bearing ringwoodite (b,d). First principles calculations 
within LDA (solid lines) are compared to available experimental data (symbols) at P=0 
GPa. 
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by these two factors. 

The temperature dependence of elastic moduli and sound velocities of 
Fe-free and Fe-bearing ringwoodite at ambient pressure are shown in Fig . 



2] Fo r x — 0, the agreement between experimental results (IJackson et al 



20001 ) and our findings is outstanding (Figs. [2h and [2b). Similarly, our 



predicted aggregate properties for x = 0.125 are in excellent correspondence 



142 with experiments having x = 0.09 (IMayama et al. 



2005 



Sinogeikin et al. 



20031 ) (Figs. [2b and[2H). Predicted Fe-bearing K is larger than Fe-free K at 



300 K (Table [T]), and dK/dT is more negative for x = 0.125 than for x = 0. 
On the other hand, G decreases with iron, but likewise K, dG/dT is more 
negative for x = 0.125 than for x = 0. Predicted compressional, shear, and 



bu 

by 



k velocities as function of temperature are smaller than those reported 



Sinogeikin et al.l (120031 ). which can be understood given the difference in 



iron content (Fig. [2H). The x dependence of elastic moduli and velocities 
for a—, (3—, and 7— (Fe^Mg^^SiC^ at high temperatures and pressures 
is shown in Fig. [31 We find dK/dx to be positive for all three phases, 
while dG/dx, dVp/dx, dVs/dx, and dV^/dx are negative. For small x, a 
linear trend given by our results of elastic moduli and velocities compares 
well to experimental data of olivine and ringwoodite. On the other hand, 
experimental values of wadsleyite are more scattered and deviate the most 
from the proposed linear behavior. Detailed dependence on pressure and 
temperature of dK/dx, dG/dx, dVp/dx, dVs/dx, and dV^/dx are shown in 
Figs. [Hand [51 dK/dx and dG/dx for all three phases exhibit qualitatively 
similar behavior in the pressure range considered. At low pressure they are 
quite sensitive to temperatures, but they seem to converge at high pressure 
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X 

Figure 3: (Color online) Dependence on low iron content at P = GPa and T = 300 K of 



a) elastic modul i, K and G, and b) velo cities, Vp and V s, (lines) compared to ex pe rimental 
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Figure 4: (Color online) Pressure dependence of dK/dx and dG/dx for a) olivine, b) 
wadsleyite, and c) ringwoodite. 
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Figure 5: (Color online) Pressure dependence of dV/dx for a) olivine, b) wadsleyite, and 
c) ringwoodite. 



161 (Figs. @Ji and St). dVp/dx, dVs/dx, and dV$/dx for all three phases are also 

162 more sensitive to temperature at lower pressures. 



163 4. Geophysical implications: The 520 km Discontinuity 

164 The seismic discontinuity near 520 km depth is often attribut e d to the 



165 phase change of wadsleyite to ringw oodite (IKatsura and 



1990 



3, 



1989 



Shearer! . 



Revenaugh and Jordan! . 1 199 ll ). It is more likely to be unobserved than 



167 either of its near neighbors at 4 10-km and 660-km depth. It i s, on average, a 

168 smaller amplitude feature (e.g. iRevenaugh and Jordan! . Il99ll ) such that less 

169 frequent observation is to be expected. In som e studies the 520 km disconti 



nuity appears as a spli t arrival or doublet (e.g 



Deuss and Woodhouse 



2001 



Chambers et al 



2005 



Bagley et al 



20091 ). When split, the two discontinu- 



172 ities are observed at depths of approximately 500 and 560 km (IDeuss and Woodhouse 



200ll ). Notably the sum of the two seismic features is larger than typical 
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174 non-split observations. Whether this is the result of an upward bias in iden- 

175 tifying split arrivals or the result of greater net velocity contrast it is not 

176 clear. Therefore accurate data on elasticity of wadsleyite and ringwoodite 

177 are critical for investigating the role of the transformation /3 — > 7 on the 520 

178 km seismic discontinuity. W e use our results on aggregate properties of (3— 



179 ( Nunez Valdez et al. 



20121 ) and 7— (Fea;,Mgi_ x .)2Si04 at temperatures and 

wo pressures encompassing the TZ to estimate the magnitude of the discontinu- 

181 ity across the phase transition. Although this is a divariant phase transition 

182 and calculation of the two-phase loop is beyond the scope of this work, we 

183 can clearly calculate velocity increases throughout the entire transition. As 

184 we saw in the previous section, experimental studies dealing with simultane- 

185 ous high pressures and temperatures offer limited data, and analyses usually 
we extrapolate results at ambient conditions either in temperature or pressure 

187 to conditions near 520 km depth (~18 GPa and ~1600 K). The lack of other 

188 source of knowledge makes it difficult to outline conclusions and/or explain 

189 the nature of the 520 km seismic discontinuity. With this paper we hope to 

190 advance the understanding of this discontinuity. 

191 To quantify the magnitude of the discontinuity across the to 7 transition 

192 at finite temperatures we use the contrast A of a particular property M 

193 defined as: 

AM = = Mx f ] x 100, (15) 

2 

194 where M could be density, elastic modulus, or velocity. Table [2] and Fig. 

195 [6] show our calculated contrasts at finite temperatures for the (3 — » 7 tran- 

196 sition in Fe-free and Fe-bearing phases. We first notice that Ap is almost 

197 independent of temperature and pressure with iron having an insignificant 

15 



(3^y Mg 2 Si0 4 p^Y (Fe x ,Mg 1x ) 2 Si0 4 



o 

5 
15 

°r 10 

"55 

re 







(a.1) — 300 K 

1000 K 

1500 K 

1700 K 

1673 K R1992 
+ 1500 K Y2008 

AG+5 ° SL2006 


(b.1) x=0.125: 300 K 

1000 K 

1500 K 

1700 K 

^ 1673 K, x=0.12 

AG+5 


. A K 


AK 


^ Ap 
P 


Ap 










AV +10 ^ 

_ A V 5 


A V +10 


IE 
o 

. AV D 


AV s +5 


(a.2) 5 


AV 

(b.2) 



16 17 18 19 20 21 16 17 18 19 20 21 
Pressure (GPa) Pressure (GPa) 



Figure 6: (Color online) Density, elastic, and velocity contrasts (lines) compared to lab- 
oratory and seismic d a ta acr oss the Fe -free (a,c) andFe-bearing ( b,d) — > 7 transition. 



R1992: 



Rigden et al 



(120061) : x = 0.12: 



(1992); Y2008: 



Sinogeikin et al 



Yu et al 



(120031) 



(120081) : SL2006: 



Lawrence and Shearer 
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effect (Figs. Ek and Eb), as suggested by lYu et al.l ( 120081 ) . Extrapolated 
experimental results and seismi c data also suggested that iron should have 



200 a small effect on this quanti ty (IRigden et al. 



2006 



Sinogeikin et al 



1991 



Lawrence and Sheared . 



20031 ). For x = we o btain a very good agreement 



with a previous estimation by 



slightly smaller than that by 



Yu et al. 



(2008) 



Sinogeikin et al. 



For x = 0.125 our result is 



(120031 ) for x = 0.09. Contrasts 



of elastic moduli (Figs. Ek and [6b) indicate that AG is more sensitive to 
temperature than AK. This dependence is greater in x = 0.125 than in 
x = 0. In the 16-21GPa range, both AG and AK slightly decrease with 
pressure, irrespective of x. Figs. Efc and Eli show that velocity contrasts, 
AVp, AVs, and AV$, exhibit weak pressure dependence, decreasing only 
slightly with pressure. Temperature affects AVs rnore than AVp or AV$. 
For x = 0, a compar i son of our results to extrapolations of experimental 



2ii data by 



Rigden et al. 



(1991) to 18 GPa and 1673 K and seismic data by 



Lawrence and Shearerl (120061 ) shows predicted contrasts to b e larger. AV$ is 



213 in better agreement with a previous pre diction at 1700 K (lYu et al. 



20081 ) 



and with a contrast from seismic studies (lLawrence and Shearerl . 120061 ) . Ve- 

215 locity contrasts increase with increasing x, however AVs is the most enhanced 

216 by x, followed by AVp, and AV$. Extra polations of AVp and AV s to 18 GPa 

217 and 1673 K with x = 0.12 reported by ISinogeikin et al.l (120031 ) are smaller 

218 than our predicted contrasts, though AVs is closer to our prediction. 

219 In a n attempt to explai n the presence of two mid-transition zone disconti- 



220 nuities 



Saikia et al 



(120081 ) investigated the solubility of CaSiOs perovskite in 
majorite garnet at high-pressure (15-24 GPa) and high-temperature (1400° 
C-1600 C). They concluded that in fertile peridotite (i.e., peridotite en- 
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/MMgi-^Fe^SiC^ 



V (A 3 ) 


K (GPa) 


G (GPa) 


X 


Reference 


541.3 


164.4 


107.7 





Nunez Valdez et al. (20121, Single crvstal/DFT 




170 


108 





Li et al. ("19961. Polv-crvstal /US 


535.8(0.2) 


170(2) 


115(2) 





Zha et al. (19971. Single crystal /BS 




169.7(1.9)-170.7(2) 


113.9(0.7)-114.1(0.8) 





Isaak et al. (20071, Polv-crvstal /RUS 


541.5 


169.5 


101.7 


0.125 


This study, Single crystal/DFT 


539.4(4) 


170(3) 


108(2) 


0.08 


Sinogeikin et al. (19981, Single crvstal/BS 




170.8(1.2) 


108.9(0.4) 


0.08 


Isaak et al. (20101, Polv-crvstal /RUS 




165.72(6) 


105.43(2) 


0.09 


Mavama et al. (20041. Polv-crvstal /RUS 




172(2) 


106(1) 


0.12 


Li and Licbcrmann (2000), Poly-crystal/US 




171.3(3) 


108.7(2) 


0.13 


Liu et al. (20091. Polv-crvstal/US 


7-(Mgi_ a; ,Fo I )2Si04 


V (A 3 ) 


K (GPa) 


G (GPa) 


X 


Reference 


527.6 


184.4 


120.95 





This study, Single crystal/DFT 


525.3 


184 


119 





Weidner et al. (19841. Single crvstal/BS 


525.3 


185(3) 


120.4(2) 





Jackson et al. (20001, Single crvstal/BS 




185(2) 


120(1) 





Li (20031, Poly-crystal/US 




185(2) 


127(1) 





Higo et al. (20061, Polv-crvstal/US 


527.7 


186.3 


115.3 


0.125 


This study, Single crystal/DFT 


526.2(4) 


188.3(30) 


119.6(20) 


0.09 


Sinogeikin et al. (20031, Single crvstal/BS 




185. 11(0.16)-185. 17(0.17) 


118.27(0.06) 


0.09 


Mavama et al. (20051. Polv-crvstal /RUS 




187(2) 


116(1) 


0.20 


Higo et al. (20061. Polv-crvstal/US 



Table 1: Results on wadsleyite and ringwoodite for volume (V), bulk (K), and shear (G), 
moduli at ambient pressure and temperature. US: Ultrasonic techniques; BS: Brillouin 
scattering techniques; RUS: Resonant Ultrasonic techniques. 
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riched in Ca and Al) at 1400° C, the wadsleyite to ringwoodite phase change 
produces a strong discontinuity at ~500-520 km depth, while the exsolution 
of Ca-perovskite produces a weak discontinuity near 540 km. At 1600° C, the 
two merge to form a single discontinuity at 540-560 km. In MORB-like com- 
positions, wadsleyite and ringwoodite are effectively absent, b ut exsolution 
of Ca- perovskite causes a velocity discontinuity near 560 km. 



Saikia et al. 



(120081 ) concluded that a mechanical mixture (or seismically averaged assem- 
blage) of peridotite and MORB would have two discontinuities: one near 500 
km due to the wadsleyite to ringwoodite phase change, and a second near 
560 km due to the exsolution of Ca-perovski te from garnet, consistent with 



233 a hypothesis of 



Deuss and Woodhousd (120011 ). 



Of some question here is the association of the MORB-like mantle with 
the non-MORB component. If MORB is associated with peridotitic man- 
tle, the strength of the 500 km discontinuity would decrease with respect to 
that of the 520 km proportionally to the percentage of MORB. Since AVs is 
smaller for the ex-solution of C a-perovskite from m ajorite garnet in MORB- 



239 like material by a factor of 2 (jSaikia et al 



20081 ). the summed strength of 



the two seismic discontinuities, at 500 km and 560 km, in a MORB plus 
peridotite mixture would b e smaller than that of the 52 km alone. This 



is opposite to observations (IDeuss and Woodhouse 



20011 ) . and whether this 



is the result of an upward bias in identifying split arrivals or the result of 
greater net velocity contrast is not clear. If MORB is associated with MORB- 
depleted mantle (harzburgite or dunite) the greater wadsleyite fraction in the 
depleted mantle might compensate and increase in AVs for the 500 km dis- 
continuity. However, the smaller fraction of iron in MORB-depleted mantle 
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T(K) 


1500 


1700 


X 


0.0 


0.125 


0.0 


0.125 


Ap 


1.89 


1.92 


1.91 


1.94 


AK 


7.94 


8.31 


8.13 


8.41 


AG 


8.14 


10.26 


8.42 


10.55 


AV P 


3.07 


3.58 


3.17 


3.67 


AV S 


3.12 


4.17 


3.26 


4.31 


AV& 


3.03 


3.20 


3.11 


3.23 


A(pV P ) 


4.96 


4.99 


5.08 


5.11 




5.02 


5.04 


5.17 


5.20 



Table 2: Predicted contrasts in % across the j3 — > 7-(Mgi_ a; ,Fe a ;)2Si04 transition at 18 
GPa. 



248 fjJaques and Green 



19801 ) should decrease AVs (as shown in Table [2]) of the 

249 to 7 transition. The extent to which these effects counter-act each other 

250 and whether the greater summed strength of the two discontinuities can be 

251 explained by this picture, remains an open question. 



252 5. Conclusions 

253 For the first time, we have presented parameter-free first-principles results 

254 of high pressure and high temperature aggregate elastic properties and sound 

255 velocities of Fe-bearing ringwoodite. We used the QHA and a novel metho d 



256 of calculating elasticity at high temperatures (iWu and Wentzcovitchl . 120 111 ). 

257 Treatment of strain Griineisen parameters via isotropic averages reduced 

258 greatly the computational cost of the task, which, otherwise, would have been 
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much more intensive and lengthy. Within the QHA limit, our predictions for 
elastic and acoustic properties for x = were found to be in very good 
agreement with available exper i mental data at 30 K and ambient pressure 



Li et al 



1996 



Weidner et al 



Zha et al. 



1997 



Isaak et al. 



20071 : 



Li 



2003 



1984 ) and molecular dynamics simulations (ILi et al 



Higo et al 



2006; 



20061 ) at 



264 high pressures and temperatures. For x = 0.125 our resul ts compared very 



265 well with experimental data in the range x = 0.08 — 0.2 ( Sinogeiki n et al 



1998 



Li and Liebermann 



2000 



Liu et al. 



2009 



Sinogeikin et al. 



2003 



Higo et al. 



20061 ). High-temperature and ambient pr essure resu 



2004 



Isaak et al. 



20K 



Jackson et al 



repro duced experiment al trends very well (llsaak et al 



2000 



t s in ringwoodite also 



2007 



Mayama et al 



Mavama et al 



2005 



Sinogeikin et al. 



20031 ) for Fe-free and Fe-bearing samples. 

Overall our predictions showed well defined changes in the elastic and 
acoustic properties of the /3— and 7— phases near conditions of the 520 km 
seismic discontinuity. We show that pressure tends to decrease contrasts 
across the j3 — > 7 transition while temperature and iron concentration tend 
to enhance them. The absence of global observations of the 520 km discon- 
tinuity could suggest regions of the TZ with less iron and/or smaller olivine 
content and irregular temperature topography. Other considerations to try 
to explain the intermittent nature of the discontinuity would involve changes 
in the pyroxene / garnet / Ca-pv system and the amount of water present in the 
TZ. These issues will be addressed in future similar studies including other 
relevant phases. 
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